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Abstract 

The extent of a locally equilibrated parton plasma in d -|- An collisions at yCMv = 200 
GeV is investigated as a function of collision centrality in a nonequilibrium-statistical 
framework. Based on a three-sources model, analytical solutions of a relativistic diffusion 
equation are in precise agreement with recent data for charged-particle pseudorapidity 
distributions. The moving midrapidity source indicates the size of the local thermal 
equilibrium region after hadronization. In central d -|- An collisions it contains about 
19% of the produced particles, and its relative importance rises with decreasing centrality. 
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The production and identification of a transient quark-gluon plasma in local thermal 
equilibrium is of basic importance in relativistic heavy-ion physics. In this Letter we 
propose nonequilibrium-statistical methods to investigate analytically the gradual ther- 
malization occuring in the course of particle production at the highest available energies. 
The approach is tailored to identify the fraction of produced particles in local thermal 
equilibrium from their distribution functions in pseudorapidity. It yields indirect evidence 
for the extent of a locally equilibrated parton plasma. 

Recently pseudorapidity distributions of primary charged particles have become avail¬ 
able PP as functions of centrality in d -|- Au collisions at a nucleon-nucleon center-of-mass 
energy of 200 GeV. They are investigated within a nonequilibrium-statistical framework 
that is based on analytical solutions of a Relativistic Diffusion Model (RDM). 

Our investigation is based on a linear Fokker-Planck equation (FPE) for three com¬ 
ponents Rk{y,t) of the distribution function in rapidity space |21E1I11 
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with the rapidity j/ = 0.5 ■ ln{{E + p)/{E — p)). The diagonal components Dy of the 
diffusion tensor contain the microscopic physics in the respective Au-like (k=l), d-like 
(k=2) and central (k=3) regions. They account for the broadening of the distribution 
functions through interactions and particle creations. In the present investigation the 
off-diagonal terms of the diffusion tensor are assumed to be zero. The rapidity relaxation 
time Ty determines the speed of the statistical equilibration in y-space. 

As time goes to infinity, the mean values of the solutions of Eqs. m approach the 
equilibrium value yeq- We determine it from energy- and momentum conservation jSl El 
in the system of Au- and d-participants and hence, it depends on impact parameter. 
This dependence is decisive for a detailed description of the measured charged-particle 
distributions in asymmetric systems: 

_ . /o / < > ^xp{ymax)+ < ml{h) > exp{-ymax) 

^ ^ m^{b) > exp{y^ax)+ < > exp{-y^a,) ^ ^ 


with the beam rapidities jb = ±2/maa:, the transverse masses < >= 

2 (^)+ < Pt >^), and masses mi^ 2 (&) of the Au- and d-like participants that depend 
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on the impact parameter b. The average numbers of participants iVi 2 (&) in the incident 
gold and deuteron nuclei are calculated from the geometrical overlap. The results are 
consistent with the Glauber calculations reported in P which we use in the further anal¬ 
ysis. The corresponding equilibrium values of the rapidity vary from j^q = - 0.169 for 
peripheral (80-100%) to j^q = - 0.944 for central (0-20%) collisions. They are negative due 
to the net longitudinal momentum of the participants in the laboratory frame, and their 
absolute magnitudes decrease with impact parameter since the number of participants 
decreases for more peripheral collisions. 

The RDM describes the drift towards i/eq in a statistical sense. Whether the mean val¬ 
ues of the distribution functions Ri and i ?2 actually attain j/eg depends on the interaction 
time Tint (the time the system interacts strongly, or the integration time of (P)). It can be 
determined from dynamical models or from parametrizations of two-particle correlation 
measurements. For central An -|- An at 200 A GeV, this yields about Tint — 10/m/c P, 
which is too short for Ri and R 2 to reach equilibrium. Note, however, that this does not 
apply to Req which is born near local equilibrium at short times (in the present calcula¬ 
tion, at t = 0 due to the 5—function intitial conditions), and then spreads in time through 
diffusive interactions with other particles at nearly the same rapidity. 

Our analytical diffusion model is consistent with, and complementary to parton cas¬ 
cade models where stopping involves large sudden jumps in rapidity from hard scatterings 
(eg- 0 ), because even hard partons can participate signihcantly in equilibration processes, 
as is evidenced by the high-pr suppression found in An -|- An at RHIG. 

Nonlinear effects are not considered here. Their possible role in the context of rela¬ 
tivistic heavy-ion collisions has been discussed in pcnip. These account to some extent 
for the collective expansion of the system in y—space, which is not included a priori in 
a statistical treatment. In the linear model, the expansion is treated through effective 
diffusion coefficients DyR that are larger than the theoretical values calculated from the 
dissipation-fluctuation theorem that normally relates Dy and Ty to each other m- One 
can then deduce the collective expansion velocities from a comparison between data and 
theoretical result. 

The FPE can be solved analytically in the linear case with constant Dy. For net- 
baryon rapidity distributions, the initial conditions are (5-functions at the beam rapidities 
Ub = ±ymax- However, it has been shown that in addition there exists a central (k=3. 
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equilibrium) source at RHIC energies which accounts for about 14% of the net-proton 
yield in Au -I- Au collisions at 200 AGeV |3], and is most likely related to deconhnement. 
For d -|- Au, net-proton rapidity distributions are not yet available. 

For produced particles, the initial conditions are not uniquely defined. Our previ¬ 
ous experience with the Au -|- Au system regarding both net baryons jl], and produced 
hadrons [121 favors a three-sources approach, with ^-function initial conditions at the 
beam rapidities, supplemented by a source centered at the equilibrium value yeq. This 
value is equal to zero for symmetric systems, but for the asymmetric d -|- Au case its 
deviation from zero according to m is decisive in the description of particle production. 

Physically, the particles in this source are expected to be generated mostly from gluon- 
gluon collisions since only few valence quarks are present in the midrapidity region at 
■\/snn — 200 GeV |^. Particle creation from a gluon-dominated source, in addition to the 
sources related to the valence part of the nucleons, has also been proposed by Bialas and 
Gzyz PI- The final width of this source corresponds to the local equilibrium temperature 
of the system which may approximately be obtained from analyses of particle abundance 
ratios, plus the broadening due to the collective expansion of the system. Formally, the 
local equilibrium distribution is a solution of m with diffusion coefficient = Dy'^, and 
(5-function initial condition at the equilibrium value. 

The PHOBOS-collaboration has analyzed their minimum-bias data successfully using 
a triple Gaussian fit M- This is consistent with our analytical three-sources approach, 
although additional contributions to particle production have been proposed PI- Beyond 
the precise representation of the data, however, the Relativistic Diffusion Model offers an 
analytical description of the statistical equilibration during the collision and in particular, 
of the extent of the moving midrapidity source which is indicative of a locally equilibrated 
parton plasma prior to hadronization. 

With (5—function initial conditions for the Au-like source (1), the d-like source (2) 
and the equilibrium source (eq), we obtain exact analytical diffusion-model solutions as 
an incoherent superposition of the distribution functions Rk{y,t) because the differential 
equation is linear. The three individual distributions are Gaussians with mean values 

< yi,2{^) >= 2/eg[l - exp{-t/ry)] =F ymax exp {-t/Ty). (3) 

for the sources (1) and (2), and yeq for the moving equilibrium source. Hence, all three 
mean values attain yeq{b) as determined from (|2)) for t—> oo, whereas for short times the 
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mean rapidities are smaller than, but close to the An- and d-like values in the sources 1 
and 2. The variances are 


a 


1,2,eg 


it) = -exp(-2t/rj]. 


(4) 


The charged-particle distribution in rapidity space is then obtained as incoherent 
superposition of nonequilibrium and local equilibrium solutions of o 


T.„,) ^ + N;My,T,„) + NSIi‘:^(v,T,„,) (5) 

with the interaction time Tint (total integration time of the differential equation). In the 
present work, the integration is stopped at the value of Tint/Ty that produces the minimum 
with respect to the data and hence, the explicit value of Tint is not needed as an input. 
The result for central collisions is Tint/Ty ~ 0.4. 

The average numbers of charged particles in the Au- and d-like regions are pro¬ 
portional to the respective numbers of participants A^i, 2 , 


N, 


1,2 
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with the constraint Ad))* = Ad^ -|- Ad^j -|- Here the total number of charged particles 
in each centrality bin A"*))* is determined from the data. The average number of charged 
particles in the equilibrium source is a free parameter that is optimized together with 
the variances and TintjTy in a y^-ht of the data using the CERN minuit-code jinj. With 
known including its dependence on centrality, one could then determine Ty and Dy, 
but this is beyond the scope of the present work. 

The result of the RDM calculation is shown in Fig. 1 for central collisions (0-20%) of 
d -|- Au. The charged-particle yield is dominated by hadrons produced from the Au-like 
source, but there is a sizeable equilibrium source that is more important than the d-like 
contribution. This thermalized source is moving since y^q has a negative value for d -|- 
Au, whereas it is zero for symmetric systems. 

The total yield is compared to PHOBOS data [T] which refer to the pseudorapidity 
1 ] = —ln[tan{9/2)] since particle identihcation was not available. As a consequence, there 
is a small difference to the model result in y-space {y ~ rj) which is most pronounced in 
the midrapidity region. It is removed when the theoretical result is converted to r^-space 
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through the Jacobian 

J{rj, {m)/{pT)) = cosh(? 7 ) ■ [1 + {{m)/{ pt)Y + sinh^(? 7 )]"^/^. (7) 

Here we approximate the average mass < m > of produced charged hadrons in the central 
region by the pion mass and use a mean transverse momentum < pt > = 0.4 GeV/c. 
In the Au-like region, the average mass is larger due to the participant protons, but 
since their number Zi < 5.41 is small compared to the number of produced charged 
hadrons in the d + An system, the increase above the pion mass remains small: < m >~ 
nip-Zi/— Zi)/Nl^ ^ O.lTGeH. This increase turns out to have a negligeable 
effect on the results of the numerical optimization, where we use < m > / < Pt >= 0.45 
for the Jacobian transformations in the three regions. For reasonable deviations of the 
mean transverse momentum from 0.4 GeV/c, the results remain consistent with the data 
within the experimental error bars. 

The equilibrium source in the light and asymmetric d + An system is found to contain 
only 19% of the produced charged hadrons in central collisions. A previous result ^21 
for An + An in the three-sources-RDM shows that the equilibrium source for particle 
production tends to be larger in the heavy system 1121 at the same energy. Note, however, 
that the results of the minimization are not unique for An + An due to the symmetry 
of the system. 

The results for the mean values and variances of produced charged hadrons in d 
+ An collisions as functions of impact parameter are shown in Fig. 2. Here the 
average impact parameters for the five centrality cuts k are determined according to 
< bk >= J bak{b)db/ J ak{b)db with the geometrical cross sections (Jk{b) in each bin. 
Whereas the total particle number and the particles created from the Au-like source de¬ 
crease almost linearly with increasing impact parameter, the magnitude of the equilibrium 
source is roughly independent of centrality. As a consequence, particle production in the 
equilibrium source is relatively more important in peripheral collisions. The variance of 
the central source lies for sufficiently small impact parameters between the values for the 
An- and d-like sources. 

The model calculations are converted to r^-space and compared with PHOBOS data 
for hve centrality cuts PP and minimum bias jT3] in Fig. 3. The minimization procedure 
yields precise results so that reliable values for the relative importance of the three sources 
for particle production can be determined. Table 1. The rapidity relaxation times and 
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diffusion coefficients can also be obtained from 0,0, but this requires an independent 
information about the interaction times. A small discrepancy in case of the most pe¬ 
ripheral collisions (80-100%) is a consequence of the three straggling data points in the 
region -4 < r; < — 3. The observed shift of the distributions towards the Au-like region 
in more central collisions, and the steeper slope in the deuteron direction as compared 
to the gold direction appear in the Relativistic Diffusion Model as a consequence of the 
gradual approach to equilibrium. 

Given the structure of the underlying differential equation that we use to model the 
equilibration, together with the initial conditions and the constraints imposed by Eqs. 
o and 0, there is no room for substantial modihcations of this result. In particular, 
changes in the impact-parameter dependence of the mean values in 0 that are not in 
accordance with 0 vitiate the precise agreement with the data. 

To conclude, we have investigated charged-particle production in d -|- An collisions 
at ■\/sisiN= 200 GeV as function of centrality within the framework of an analytically 
soluble three-sources model. Excellent agreement with recent PHOBOS pseudorapidity 
distributions has been obtained, and from a x^-minimization we have determined the 
diffusion-model parameters very accurately. 

For central d -|- An collisions, a fraction of only 19% of the produced particles arises 
from the locally equilibrated midrapidity source. Although this fraction increases to¬ 
wards more peripheral collisions, the formation of a thermalized parton plasma prior to 
hadronization can probably only be expected for more central collisions. 

The d -|- An results show clearly that only the midrapidity part of the distribution 
function reaches thermal equilibrium, whereas the interaction time is too short for the 
d- and Au-like parts to attain the thermal limit. The same is true for the heavy An -|- 
Au system at the same energy, but there the precise fraction of particles produced in the 
equilibrium source is more difficult to determine due to the symmetry of the problem. 

One of the authors (GW) acknowledges the hospitality of the Faculty of Sciences at 
Shinshu University, and financial support by the Japan Society for the Promotion of 
Science (JSPS). 
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Table 1. Produced charged hadrons as functions of centrality in d + Au collisions at 
= 200 GeV, Yb = ± 5.36 in the Relativistic Diffusion Model. The variance of the 
central source in y—space is The number of produced charged particles is for the 
sources 1 and 2 and for the equilibrium source, the percentage of charged particles 
produced in the thermalized source is 


Centrality{%) 

< 

Nk 

N!, 

NS 

n:u%) 

0-20 

3.99 

131 

19 

35 

19 

20-40 

3.95 

78 

17 

31 

25 

40-60 

5.70 

33 

11 

38 

46 

60-80 

7.44 

9 

5 

35 

71 

80-100 

6.89 

2 

2 

24 

86 

min. bias 

4.04 

56 

15 

21 

23 
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Figure captions 


Fig. 1. Charged-particle rapidity spectrum dNch/dy in the Relativistic Diffusion Model 
(RDM) for central (0-20%) d -|- An collisions at ^/s]^ = 200 GeV, upper curve. 
Dashed curves are contributions from the Au- and d-like sources, respectively. The 
dashed area is the yield from the moving thermalized central source of partonic 
origin, cf. text. PHOBOS pseudorapidity data dNch/dr] (r; ~ ?/) P are shown to 
illustrate the difference in y— and 77 —space. See Fig. 3 for a fit with the proper 
Jacobian transformation. 

Fig. 2. Mean values (upper part) and variances (lower part) of the three sources for pro¬ 
duced charged particles in rapidity space as functions of the mean impact parameter 
for five centrality bins, symbols. Lines are drawn to guide the eye. Dashed lines 
have Au- and d-like sources, solid lines correspond to the thermalized central source. 
The top solid line gives the mean total number of produced charged particles. 

Fig. 3. Calculated pseudorapidity distributions of charged particles from d -|- Au colli¬ 
sions at y'SNN = 200 GeV for five different collision centralities, and minimum-bias 
in comparison with PHOBOS data Pd. The analytical RDM-solutions are op¬ 
timized in a fit to the data. The corresponding minimum y;^-values (top left to 
bottom right) are 4.7, 5.9, 2.4, 1.7, 1.9, 2.1. See Fig. 2 and text for the resulting 
diffusion-model parameters. 
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